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Abstract 

The bound state (energy spectrum and two-spinor wave functions) solutions of the Dirac equation 
with the Hulthen potential for all angular momenta based on the spin and pseudospin symmetry are 
obtained. The parametric generalization of the Nikiforov-Uvarov method is used in the calculations. 
The orbital dependency (spin-orbit and pseudospin-orbit dependent coupling too singular 1/r^) of 
the Dirac equation are included to the solution by introducing a more accurate approximation 
scheme to deal with the centrifugal (pseudo-centrifugal) term. The approximation is also made 
for the less singular 1/r orbital term in the Dirac equation for a wider energy spectrum. The 
nonrelativistic limits are also obtained on mapping of parameters. 
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I. INTRODUCTION 



The spin or pseudospin symmetry [1,2] investigated by the framework of the Dirac equa- 
tion is one of the most interesting phenomena in the relativistic quantum mechanics to 
explain different aspects for nucleon spectrum in nuclei. This is mainly studied for the ex- 
istence of identical bands in superdeformed nuclei in the framework of a Dirac hamiltonian 
with attractive scalar S{r) and repulsive vector V{r) potentials [3]. The pseudospin sym- 
metry is based on the small energy difference between single-nucleon doublets with different 
quantum numbers and the Hamiltonian of nucleons moving in the relativistic mean field 
produced by the interactions between nucleons. The relativistic dynamics are described by 
using the Dirac equation only [4]. 

Ginocchio [5] found that the pseudospin symmetry concept in nuclei occurs when 5'(r) 
and V{r) potentials are nearly equal to each other in magnitude but opposite in sign, i.e., 
S{f) ~ —V{r) and hence their sum is a constant, i.e., E(r) = V{r) + S{r) = Cps- A 
necessary condition for occurrence of the pseudospin symmetry in nuclei is to consider the 
case E(r) = [5-7]. Further, Meng et al [8] showed that the pseudospin symmetry is exact 
under the condition of dTi{r)/dr — 0. Lisboa et al [9] studied the generalized harmonic 
oscillator for spin-1/2 particles under the condition E(r) = or A(r) = V{f) — S{r) = 0. 
The Dirac equation has been solved numerically [10,11] and analytically [4,12,13] for nucleons 
that are moving independently in the relativistic mean field in the presence of the pseudospin 
symmetric scalar and vector potentials. The exact analytical solutions of the Dirac equation 
gives the bound-state energy spectra and spinor wave functions [14,15]. 

The aim of this paper is to present an analytical bound state solutions of the Dirac 
equation for the Hulthen potential in the presence of the exact pseudospin (spin) sym- 
metry using a new approximation scheme to deal with the pseudo-centrifugal (centrifugal) 
potential term for / > (/ > 0) case. To obtain a general solution for all values of the pseu- 
dospin (spin) quantum numbers, the pseudospin (spin) symmetry and orbital dependency, 
pseudospin-orbit (spin-orbit) dependent coupling are included to the lower component of 
the Dirac equation as an integer quantum number. This component has the structure of 
the Schrodinger-like equations with the pseudo-centrifugal (spin-centrifugal) kinetic energy 
term and its solution is analyzed by using some algebraic methods and effective approaches. 
For small values of the radial coordinate r, this effective potential gives a centrifugal energy 
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term in the first approximation. The Dirac equation for the Hulthen potential is arranged 
under the condition of the exact pseudospin (spin) symmetry and it's solution is obtained 
systematically by using the Nikiforov-Uvarov (NU) method [16]. As an apphcation of the 
Dirac-Hulthen problem with the pseudospin (spin) symmetry, the relativistic eigenvalue 
spectrum for various degenerate states is presented for several pseudo-orbital (spin-orbital) 
and pseudospin (spin) quantum numbers. 

The paper is organized as follows. In Sec. 2, the outline of the parametric generalization 
of the NU method is presented. Section 3 is devoted for the relativistic quantum mechanics 
(Dirac and Klein-Gordon equations) and the additional coupling of the space scalar and vec- 
tor potentials to free particle wave equations. In Sec. 4, the pseudospin and spin symmetry 
Dirac equation in 3 -|- 1 dimensions with coupling is solved for the Hulthen potential 
using an improved approximation scheme to deal with the too singular pseudo-centrifugal 
(centrifugal) kinetic energy term «;(«;± l)/r^. The parameteric generahzation of the NU 
method is followed to obtain the energy eigenvalues and the corresponding two-spinor wave 
functions. In Sec. 5, we solve the Dirac equation with an accurate proper approximation 
made for the less singular coupling 1/r vector potential to extend the validity of the results 
to a wider range energy spectrum. Results and conclusions are presented in Sec. 6. 



II. METHOD OF ANALYSIS 



The Schrodinger-like equation including the centrifugal barrier and/or the spin-orbit cou- 
pling term has not been solved exactly for the exponential-type potentials such as Morse, 
Hulthen, Woods-Saxon, etc [4,12-15,17]. The exact solution of the Schrodinger equation 
for the exponential-type potentials has been obtained for / = 0, however, any /-state so- 
lutions have been given approximately by using some analytical methods under a certain 
number of restrictions [4,18-20]. One of the calculational tools utilized in these studies is 
the Nikiforov-Uvarov (NU) method. This technique is based on solving the hypergeometric 
type second-order differential equations by means of the special orthogonal functions [21]. 
For a given potential, the Schrodinger or Schrodinger-like equations in spherical coordinates 
are reduced to the second-order differential equation of hypergeometric type with an appro- 
priate coordinate transformation r ^ s and then they are solved systematically to find the 
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exact or particular solutions. The NU method is briefly outlined here: 

<W + "-^f'nir) + ^fnir) - 0, (1) 

where a{r) and 5(r) arc polynomials, at most, of second- degree, and r(r) is a first-degree 
polynomial. In order to find a particular solution for Eq.(l), let us decompose the wave 
function ■0n('") ^ follows: 

V'„(r) = 0(r)y„(r), (2) 

and use 

[(j(r)p(r)]' = r{r)p{r), (3) 

to reduce Eq.(l) to the form 

a{r)y';{r) + rir)y'^{r) + Ay„(r) = 0, (4) 

with 

r(r) = r(r) + 2n{r), T'(r) < 0, (5) 

where the prime denotes the differentiation with respect to r. One is looking for a family of 
eigenvalue solutions corresponding to 

X = Xn = -nT'{r)-^n{n-l)a"{r), n = 0,l,2,---. (6) 

The |/„(r) can be expressed in terms of the Rodrigues relation: 

where Bn is the normalization constant and the weight function p{r) is the solution of the 
differential equation (3). The other part of the wave function (2) must satisfy the following 
logarithmic equation 

A^' ( r\ 'Tr(r\ 

(8) 



(f){r) a{r) 
By defining 

A; = A-7r'(r). (9) 

one obtains the polynomial 

<r) = ^ [a'{r) - ?(r)] ± ^1 [^'(r) -T{r)]' - a{r) + ka{r), (10) 



where 7r(r) is a parameter at most of order one. The expression under the square root sign in 
the above equation can be arranged as a polynomial of second order where its discriminant 
is zero. Hence, an equation for k is being obtained. After solving such an equation, the k 
values are determined through the NU method. 

We may also derive an alternative parameteric generalization from the NU method valid 
for most potential models under consideration. The first step basically begins by writting 
the hypergeometric equation [21] in general parametric form as 

[r (C3 - c^r)f iP'l^ir) + [r (cg - c^r) {c, - C2r)] iP'^{r) + {-^,r^ + i^r - i^) i;^{r) = 0, (11) 

with 

r(r) = ci - Car, a{r) = (cs - c^r) r, a{r) = -^^r^ + - Ca^ (12) 

where the coefficients Cj {i — 1,2,3,4) and the analytic expressions {j — 1,2,3) have 
to be calculated for the potential model under consideration. The second step demands 

comparing Eq.(ll) with it's counterpart Eq.(l) so that we can obtain the analytic NU poly- 
nomials, energy equation, wave functions and the relevant coefficients expressed in general 
parameteric form in Appendix A of Ref . [22] . 

III. SCALAR POTENTIAL COUPLING IN RELATIVISTIC QUANTUM ME- 
CHANICS 

The Dirac equation and the Klein-Gordon (KG) equation are wave equation mostly used 
in describing particle dynamics in relativistic quantum mechanics. These two wave equa- 
tions, for free particles, are constructed using two objects: the four- vector linear momentum 
operator = ihdn and the scalar rest mass M, allow one to introduce naturally two 
types of potential coupling. One is the gauge-invariant coupling to the four- vector potential 
{A^ {t, r )}^^Q which is introduced via the minimal substitution — > — gA^, where g is 
a real coupling parameter. The other, is an additional coupling to the space-time scalar po- 
tential S(t, 1^) which is introduced by the substitution M M + S. The term "four- vector" 
and "scalar" refers to the corresponding unitary irreducible representation of the Poincare 
space-time symmetry group (the group of rotations and translations in (3 -|- l)-dimensional 
Minkowski space-time) . Gauge invariance of the vector coupling allows for the freedom to fix 
the gauge (eliminating the non physical gauge modes) without altering the physical content 
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of the problem. Many choose to simphfy the solution of the problem by taking the space 
component of the vector potential to vanish (i.e., ^). One may write the time-component 
of the four- vector potential as qAq = V{t,l^), then it ends up with two independent poten- 
tial functions in the Dirac and KG equations. These are the "vector" potential V and the 
"scalar" potential S [23,24]. 

In the relativistic units, H — c — 1, the free Dirac and KG equations are written as 



{zYd,-M)i,j,{t,l^) = Q, 



{d^d, + M'')iljj,a{t,^) = 



(13a) 
(13b) 



respectively. The convention of summing over repeated indices is used. For particles of spin 
1/2, {7^^} are 4x4 constant matrices with the following standard representation [23]: 
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/ 
-1 



a 
-a 



where / is the 2x2 unit matrix and a are the three 2x2 hermitian Pauli spin matrices. The 

vector and scalar couplings mentioned above introduce potential interactions by mapping 
the free Dirac and KG equations as 

d 
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dt 



d 



+ • ^ - [M + S{t, -^)] I ^^(t, ^) = 0, 
'-^' + [M + 5(7^)]nV'^G(^) = 0, 



(14a) 



(14b) 



respectively. This type of coupling attracted a lot of attention in the literature due to the 
resulting simplification in the solution of the relativistic problem. The scalar-like poten- 
tial couphng is added to the scalar mass so that in case when S{r) — -\-V{r), the Dirac 
equation could always be reduced to a Schrodinger-type second order differential equation 
as we shall see in the next section. The nonrelativistic limit can be obtained by taking 
EnK. — M ~ Eni and + M ~ 2M, where \E\ <^ M. Hence, the positive energy solu- 



tion is 



E,^ 



nl 



= 0, with potential 2V, nonrelativistic energy E^i 
and tpi^) stands for either ip^i'^) or i^KoC^)- The negative energy solution correspond- 
ing to case S — —V results in a trivial non-interacting theory (free fields) with solution 
2!^^^ -I- Eni^ i^C^) — [24]. The physical meaning of introducing the scalar-like potential 
into Dirac and KG equations is making one to study the confinements in quarks when taking 



the conditions S — V, S — —V and S — rjV such that 77 7^ ±1. The last case results in 
uneven contribution of the two potentials [24]. For example, suppose S — V is a. potential 
which tends to an effective positive infinite barrier at spatial infinity for the positive-energy 
particles and will be responsible for the confinement solutions (bound states or scattering 
states). However, there is an effective infinite well for the negative- energy particles which 
cannot prevent the negative-energy particles from going to infinity [25]. This provides real 
(bound state solutions) for the positive-energy particles, but imaginary (scattering state so- 
lutions) for the negative-energy particles. Similarly, when S — V tends to negative infinity 
at spatial infinity, the positive-energy particles are not confined due to the effective poten- 
tial well [25]. Therefore, the (1 + 3)- and (1 + l)-dimensional Dirac equation with S = V, 
the confinement is impossible, i.e., there must be scattering states. If S is stronger than 
V, the confinement is permanent and if, on the contrary, the V is stronger, confinement is 
impossible due to the Klein paradox (see Ref. [25] and references therein). 



IV. DIRAC EQUATION WITH COUPLING TO r'^ SINGULAR ORBITAL TERM 



The Dirac equation of a single- nucleon of rest mass M with spherically symmetric coupling 
to an attractive scalar and repulsive vector S{f) and V{r) potentials has the following radial 
component 

'a.cP + PiMc' + S{r))] ^P^,{r) = [E^, - V{r)] ^|;^,{r), (15) 



where 



P = -ihV, a 




I 
-I 



(16) 



where P is the three momentum operators, a and (3 are the usual 4x4 Dirac matrices [23] , 
c is the speed of light in vacuum and h is the Planck's constant divided by 27r. Further, 
denotes the relativistic energy eigenvalues of the Dirac particle. For nuclei with spherical 
symmetry, S{r) = S{r) and V{f) = V{r), where r is the magnitude of f. Further, the spinor 
wave functions is 



1 



(i) 



(i) 



(17) 



where Yi{9, 0) {Yj{9, 0)) and x± the spin (pseudospin) spherical harmonic and spin wave 
function which are coupled to angular momentum j with projection m, respectively. Fn^ir) 
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and Gnnif) are the radial wave functions for the upper and lower components, respectively. 
The label k has two explanations; the aligned spin j — l-\-l/2 (si/2,P3/2, etc.) is valid for the 
case of K = — (j + 1/2) and then I — Z + 1, while the unaligned spin j = 1 — 1/2 (pi/2, 0^3/2, etc.) 
is valid for the case of k = (j ' + 1/2) and then 1 = 1 — 1. Thus, the quantum number k 
and the radial quantum number n are sufficient to label the Dirac eigenstates. The Dirac 
equation (15) can be reduced to the following two coupled ordinary differential equations 
(in the relativistic units, H = c = 1): 

(J^ + ^) Fn«(r) = {M + - A(r))G„„(r), (18) 

^ -^^ Gn.(r) = (M - K« + S(r))F„.(r), (19) 

where A(r) = V{r) — S{r) and S(r) = V{r) + S{r) are the difference and the sum potentials, 
respectively. Solving Eqs.(18) and (19) leads to a second order Schrodinger-like differential 
equation with coupling to singular term and satisfying Gn^ir), 

'^^'^ _(M + K. -A(r))(M + E(r))- ''^ ~> 




where Enn 7^ when S(r) = (exact pseudospin symmetry). Since Em^ = —M is an 
element of the negative energy spectrum of the Dirac Hamiltonian, then this relation with 
the lowe spinor component is not valid for the positive energy solution. Further, a similar 
equation satisfying ^^^(r) can be obtained as 

''^^ + '^-(M + K.-A(r))(M-K. + S(r)) ■ ^'^^^^> 




where E^k 7^ —M when A(r) = (exact spin symmetry). Since En^ = +M is an element 
of the positive energy spectrum of the Dirac Hamiltonian, then this relation with the upper 
spinor component is not valid for the negative energy solution. The exact spin symmetry 
requires {dA{r)/dr = 0, i.e., A(r) = Cg ^constant), Eq. (21) turns out to be 

d'^ k{k + 1) 



dr'^ r 



2 - (M + K« - C,)E(r) + Ei^ -M' + Cs{M- K«) j i^n«(r) = 0, (22) 

2 



where k, = I and k = —{l + l) are valid for k > and k < 0, respectively, and k{k + l)/r 
is the spin-centrifugal potential term. On the other hand, the exact pseudospin symmetry 
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requires {dE{r)/dr — 0, i.e., E(r) = Cps ^constant), Eq. (20) is reduced to the form 
(P k{k — 1 



^2 + - + C^.) A(r) + - - Cp, (M + £;„0J G'n.lr) = 0, 

(23) 

where k — l + l and k — —I are vahd for k > and k < 0, respectively, and K{K — l)/r'^ is the 
pseudo-centrifugal potential term. Therefore, the degenerate states come into existence with 
the same / but different k, generating pseudospin symmetry. The components of the wave 
function are required to satisfy the boundary conditions. That is, Fni^{r)/r — )■ (G'„k(^)A 
— )■ 0) when r — )■ oo and FnK{r)/r = {GnK{j')/r = 0) at r = hold. Note that the 
analytic solutions of the above second order differential equations require approximation to 
the orbital term k{k ± l)r~'^ that results from reduction of the original Dirac equation. For 
example, the orbital term ~ r"^ has a more singularity near r = 0. 



A. Pseudospin Symmetry Solution 



The solution of the Dirac equation (23) for the Hulthen potential demands that the 
potential A(r) is exponential in r and the pseudo-centrifugal term is quadratic in 1/r. 
Hence, the difference potential is taken as the Hulthen potential [26]: 

<5r 



A(r) = -Ao: 



:, Ao = V^o = 



where S and Aq are the screening parameters to determine the range and strength, respec- 
tively. Besides, Ze is the charge of the nuclcon [27]. This potential has been studied by 
means of the algebraic perturbation calculations based upon the dynamical group structure 
SO (2,1) [28], the NU method [29], the supersymmetry and shape invariance [30], the asymp- 
totic iteration method [31] and the Biedenharn's approach for the Dirac-Coulomb problem 
[32]. Equation (23) is analytically solvable only ior I — {k, — 1). Therefore, in order to solve 
the Dirac equation for any k or Z-state, we need to apply the following shifted approximation 
scheme near the singularity (origin) to deal with the more singular pseudo-centrifugal term, 
r-^ for the case of !> [20,29,33] 



do + 



1 



(5->0 



do + 



{5ry 



1_ {6r) 
12 



- 1 

2 



+ 



1 



{5r) 



Sr 



+ 



1)2 

{Srf 



6048 172800 



+ O {{Srf) 



(25) 
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where the approximation constraint is 5 or applying the condition (5r — > 0. It should be 
noted that the physical interpretaion in introducing the parameter do to the traditional ap- 
proximation, S'^e^^ / {e^^ — 1)^, is that when performing the power series expansion and letting 
6r — )> 0, it gives r^^ — 1/12 but not r^^ as desired. Therfore, we understand that traditional 
(conventional) approximation scheme suggested by Greene and Aldrich [34] is shifted by 
a dimensionless constant do — 1/12 from the origin. Equation (25) is the correct ansatz 
to substitute (see [33] and the references therein). The above approximation resembles 

Ki r~'^ + S^f{do), where f{do) — do — 1/12. Furtehr, it is simply the addition of traditional 
(usual) approximation plus a shifting term 5^do, i.e., fa 5 [W{r) + W'^{r)/5 + Sdo] , with 
W{r) = 6/{e^^ — 1). It is apparent from the above expansion that for small values of 5, 
the dimensionless constant d^ = 1/12. However, the approximation model used in [4] is 

— W^{r) (cf. Ref. [4] and the references therein). Figure la shows a plot of the vari- 
ation of the centrifugal orbital term with respect to r, where the screening parameter 
5 — 0.1 fm"^. We observe that the improved approximation model (solid line) works well if 
compared with r"^ (dotted-solid line). The curves in Fig. lb show that the approximation 
of is independent of the value of do. The traditional approximation is plotted as a func- 
tion for different values of do. It should be noted that old approximation in Ref. [4] has the 
ansatz ~ W'^{r) (i.e., one term function in (25)). 

Now, the substitution of Eq.(24) and Eq. (25) into Eq. (23) leads to 



GnJr) = 0, (26) 



with 



F2 , ^1- 72 ■ [^') 



Further, defining 

s = e"^'' e [0, +1], Ai=u}1-i^I + k{k-1) do, 
Bi ^2ujl-iyl + K{K- 1) {2do - 1) , el^^ujl + K{K- 1) do, (28) 
recasts Eq. (26) into the simple form 

which can be easily solved by means of the NU method or applying a short-cut procedure 
given in Appendix A of Ref. [22]. 
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The procedures begin by comparing Eq. (29) with Eq. (1) giving the polynomials: 

t{s) = 1-s, a{s) = s{l-s), a{s) = -Ais^ + Bis - el^, (30) 

and with the aid of Eqs.(12)-(14), we can obtain q = 1 (for i — 1,2, 3, 4), = Ai, ^2 = -^1 
and ^3 = e^^. In addition, the relations A1-A3 yield 

C5 = 0, C6 = --, C7 = - + Ai, cs = --Bi, cg = e„„, 

Clo — ^ , Cii — 2e„K, C12 — 2k — 1, Ci3 — e^K) C14 — (31) 

and the relations A4-A6 give the essential NU polynomials: 

7r(s) = + €nn) S, k = -uj - K {k - 1) - {2k - 1) 

t{s) = 1 + 2e„« - 2 (^enn + k+^^s, t'{s) = -2 |^e„« + k + < 0. (32) 
The eigenvalue equations (6) and (9) take the forms 

Xn^ ri^ + 2n (£„« + k) and A = -pj - ft;(2e„K + k), (33) 

respectively. In setting A = A„ or alternatively using the relation A7, we obtain the eigen- 
value equation being expressed in terms of as 



- El^ + Cps{M + K«) ^-k{k- 1) 5^do + 



where 



2 n + / + l) for K>0 



, 71 = 0,1,2,3, 



(35) 



2(n-t^ for K < 

Thus, the energy spectrum can be obtained from the following energy eigenvalue equation: 



1 + 



where 



U 



Cps + 



2AoU 



S Ni 



{Cp, + M)Ao 
5' 



SM k(k - 1 

U : ^ 

An 12 



, k{k-1)^i(^1 + l) 



0, (36) 



(37) 



The two energy solutions of the above quadratic equation are 



— 



2(A2 + iViV) 



(38) 
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For a given value of n and k (or I), the above equation provides two distinct positive and 
negative energy spectra related with or respectively. One of the distinct solutions 
is only valid to obtain the negative-energy bound states in the limit of the pseudospin 
symmetry. In the presence of exact pseudospin symmetry {Cps = 0), we finally obtain 

^ M + \N^5 ± Ni^{Nf + 1) M2 - (M + \N^5y + ^{N^ + 1)k{k - 1)5' 

■ (39) 

In this regards, states with various n and I quantum numbers having same energy spectrum 
are said to be degenerate states. 

We calculate the negative bound state energy eigenvalues [4,35,36] from Eq. (38) for 
several values of the quantum numbers n and k{1) in the pseudospin symmetry limit. They 
are displayed in Table 1. The results have been calculated by using the following choices of 
parameters: M = 5.0 fm-\ Aq = 3.40 fm'^ and Cps = -4.90 /m"^ [4]. Prom Table 1, 
one can clearly see that the degeneracy between two states in the pseudospin doublets, i.e., 
nsi/2, {n - 1)4/2 for / = 1 (/ = 0), np3/2, (n - l)/5/2 for / = 2 (/ = 1), {n - 1)5(7/2 for 

/ = 3 (/ = 2), and n/7/2, {n — l)/i9/2 for / = 4 (/ = 3), etc. Our numerical approximations 
using the new approximation scheme, Eq.(25), are compared with the ones obtained using 
Hulthen square approximation (see Eq. (34) of Ref. [4]). It is worth noting that such 
approximation schemes are usually used in literature as effective approximations to deal with 
the pseudo-centrifugal kinetic energy term in the case of i > and small r. One can easily see 
how the approximation of the energy states is sensitive and dependent on the approximation 
scheme used. Note that we have introduced a small positive shift, 5'/(/ -|- 1)/12, to the 
conventional approximation scheme [34], i.e., = 8W -|- H^', in calculating the bound 
states (real energy states and corresponding wave functions). 

Now, the corresponding wave functions calculations begin by calculating the weight func- 
tion from relation A8 as 

1 / ( t(s\ 



p{s) = ^ cxp ( / ^ds = s'^-^ (1 - sf^-' , (40) 
and the first part of the wave function: 

0(s) = exp ^ds^ = s'-- (1 - . (41) 
Further the second part of the wave function can be obtained from relation as 

Jn 

VnM = c„,s-'^-^ (1 - s)-^'^-'^ ^ [s^+'^- (1 - sT^'^-'] ~ p(f "«.2-i)(l - 2s), (42) 
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where Cuk is the normahzation constant and Pn^''^^ (x) are the Jacobi polynomials defined for 
Re(i/) > —1 and Re(//) > —1 in the interval x e [— 1,+1] .. Using GnK{s) = 4>{s)yn{s), the 
lower-spinor wave function reads 

Gn.{r) = c„« (eM-enJr)) (1 - exp(-5r))'^Pi^^"-^«-^)(l - 2exp(-5r)) 
(2en« + l)„ 



n! 



(exp(-e„«5r)) (1 - exp(-5r))'' 



with 



X 2F1 {-n, n + 2 (e„„ + /?);! + 2e„„; exp(-5r)) , > 



enK<5 = ^JM^- El^ + Cp,(E„, + M) + «:(«:- 1) ^^^o > 0. 



(43) 



(44) 



where 2F1 (— JT-, n + 2 (e„K + «);! + 2e„K; exp(— (5r)) is the hypergeometric series terminates 
for n = and thus coin orges for all \-alues of real parameters a;i > and Z > 0. When Cps — 
0, then = sJm'^ - + k{k-1) with the restriction E^^ < M'^ + k{k-1) S'^do 
is required to obtain bound state (real) solutions for both positive and negative solutions of 
En^f^ in Eq. (39). Making use of the recurrence relation of hypergeometric function 



d 
ds 



2F1 (a;6;c;s) 



(^^ ) 2F1 (a + l;6+l;c+l;s), 



(45) 



we obtain the corresponding upper component F„^re(r) from Eq. (19) as 



Funir) = br, 



K,6exp{—6r) ^ 
1 — exp(— 5r)) ""^ 



(exp(-e^^(5r)) (1 - exp{-Sr)y 
(M - + Cps) 

X 2F1 {-n, n + 2 (e„« + k);1 + ^Euk] exp(-5r)) 

n5 [n + 2 (av + e„«)] (exp(-(5r))'"'=+^ (1 - exp(-5r))'' 
(1 + 2e„«) (M - En^ + Cps) 



X 2F1 ( 1 - n, n + 2 ( e„K + K + - ) ; 2 (1 + enn) ] exp(-(5r) 



(46) 



where 6„« is the new normalization factor. Based on the exact pseudospin symmetry (i.e., 
when Cps — 0, 7^ -/W^), there are only bound negative-energy states, otherwise the upper 
spinor component Fnnir) will diverge. The energy solutions obtained from Eq.(38) for a 
given values of n and k need to be negative so that Gn^ir) and Fn^ir) are defined for the 
bound states, i.e., e„«; > 0, k > 1. 
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B. Spin Symmetry Solution 



The spin symmetry arises from S'(r) ~ ™ which the nucleon move [4]. Therefore, 
we take the sum potential equal to the Hulthen potential: 

„-6r 

E(r) = -Eo^— Eo = = ZeH, (47) 

and apply the approximation in Eq.(25) dealing with the spin-orbit centrifugal term k(k + 
l)/r^. The choice E(r) = 2V{r) — >■ V{r) enables one to restore the non-relativistic solution 
when appropriate choice of parameter transformations is being adopted [35]. Thus, Eq.(22) 
can be rewritten as 



F„,(r) = 0, (48) 



with 

(M + £;^,-C,)Eo ^_M^-El^-Cs{M-E, 
Defining the new variable and parameters. 



^2 - ?2 ' ^2 - 72 ■ i^yj 



s = e-*^ e [0, +1], A2 = ujI + uI + + do, B2 = 2cul + + k{k + 1) (2rfo - 1) , 

4k = u;^ + k (« + 1) do, (50) 

recasts Eq.(48) as 

Following the previous procedures, we obtain 

f (s) = 1 - s, (7(s) = s(l - s), ?(s) = -A2s'^ + B2S - el^, (52) 

and 

7r(s) = - (k + 1 + enn) s, k ^ - K, {k + 1) - {2k + 1) (53) 

t{s) = 1 + 2enn - 2 (^£n« + + S, (54) 

Also, the parameters A and A„ take the forms: 

Xn^n'^ + 2n (£„« + K + 1) and A = - ('^ + 1) C^^nn + K + 1), (55) 



14 



giving the energy eigenvalue equation: 



No6 



E' - 

riK 



2Eoiy 
No 



2 , TM k{k+1) 



So 
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where 



0, (56) 

(57a) 
(57b) 



{2(n + l + l) for K>0 
, n = 0,l,2,3,--- . 
2{n-l) for K < 

The two energy solutions of the quadratic equation (56) can be obtained as 

{N^Cs + 2iV2SoM/) ± 5^^{NiCs + 2iV2SoW^)' + ANi (S^ + N^S^) -™-W^^ 



2(E2 + iV|52, 



(58) 

For a given value of n and k (/), we obtain two distinct positive and negative energy spectra 
related with E^^ or E~^, respectively. However, the positive-energy solution is valid for the 
spin symmetry limit. In the presence of exact spin symmetry (C^ = 0), we can simply obtain 



E. 



± 



-M + \Nl5 ± (iV| + 1) M2 - (-M + \Nl5f + ^(iVf + 1)k{k + 1)5' 

MTi 



(59) 



Using Eq. (58), we calculate a few positive energy levels for various values of quantum 
numbers n and in the spin symmetry limit. In Table 2, wc present some numerical 
values with the following choices of parameters: M = 5.0 fm^^, Sq = 3.40 fm~^ and 
Cs = 4.90 fm^^. From Table 2, one can clearly see that the degeneracy between two states 
in the spin doublets, i.e., (npi/2, 71^3/2) for / = 1, {nd^i2,nd^i2) for I = 2, (n/5/2, n/7/2) for 
I — 3, and {ngT/2, ngg/2) for Z = 4, etc. For example, Opi/2 with n — and k = 1 (Z = 1) is 
the partner of 6^3/2 with n — and k — —2 {I — 1). 

Next, we turn into the wave functions calculations. The calculated weight function: 



Pis) 



(l-s) 



2/t+l 



(60) 



enables us to write down the second part of the wave function as 



VnM = a^.s-'^-^ (1 - s)-^'^^'^ ^ (1 - sr^'-^'] ~ p(2-^.2.+i)^i _ 2,)^ (61) 



15 



where is the normahzation constant. Furthermore, the first part of the wave function 
reads 

(f){s) = (1 - s)^+^ . (62) 
Thus, the upper component of the wave functions, F„«;(s) = becomes 



{expi-SnJr)) (1 - exp(-5r))''+^ p(2e^.M+i) (^i _ 2exp(-5r)) 



anK. ^^^""^ t^^"- {exp{-enJr)) (1 - exp{-6r))'^^^x ^Fi {-n,n + 2 {snn + k + 1) ; 1 + 2£„«; exp(-5r)) , 



where 



.5 = - E^^ + Cs{E^. - M) + « + 1) > 0. 



(63) 



(64) 



Note that the hypergeometric series 2-F1 {—n^ n + 2 (e„K + /t + 1) ; 1 + 2£:„k; exp(— (5r)) is ter- 
minated for n = and thus it converges for all values of real parameters a;2 > and k > 0. 
In case when = 0, then £„k<^ — \J ~ + k{k + 1) S^dg with a restriction for real 
bound states that < + k (k + 1) 5^do for both positive and negative solutions of 
En^ in Eq. (58). Thus, the corresponding spin-symmetric lower-component GnK{f) takes the 
form: 



(exp(-£„«,5r)) (1 - exp(-5r)) 



ft+i 



K + 1) 5exp(— 5r) ^ K 

(1 — exp(— 5r)) "'^ r 



(M + S„,-C,) 
X 2F1 (-n, n + 2 K -Fl) ; 1 2£„^; exp(-5r)) 

71(5 [n + 2 (/^ + 1 + £^^)] (exp(-^r))"-+' (1 - exp(-(5r))''+' 
(l + 2£„«) (M + Kk-C,) 



X 2F1 ( 1 - n, n + 2 ( + K + - ) ; 2 (1 + £„^) ; exp(-5r) ) , 



(65) 



where Em^ ^ —M when Cg — 0, exact spin symmetry and hnn is the normalization constant. 

Prom the above expression, we see that there are only bound positive-energy states, 
otherwise the lower spinor component Gnuij') will diverge. For a given values of n and k, we 
choose the suitable solution that makes Gnni''^) and FnK{r) satisfy the restriction conditions 
for the bound states, i.e., e„K > 0, k, > —1 and Enn are positive. 

A careful inspection of the relationship between the present set of parameters 
(a;|, A2, S2) and the previous set of parameters {uj1,i'1, Ai, Bi). provides that the spin 
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symmetric positive energy solution can be simply obtained from the pseudospin symmetric 
negative energy solution by making the replacements [37]: 

Fnnir) ^ Gnnir), V (r) -V{r) (or Sq O -Aq), k{k + 1) ^ k{k - 1) {or k± 1), 

C, O -Cp„ o -E-^, ul O u\ and vl O (66) 

That is, with the above replacements, Eqs.(38) and (43) yield Eqs.(58) and (63) and the 
vice versa is true. 

Let us now present the non-relativistic limit. This can be achieved when we set Cg — 0, 
k(k + 1) l{l-\-l),T.Q — Vq — Ze^b and using the mapping — M ~ Eni and EnK + M ~ 
2m in Eqs.(59) and (63), then energy spectrum (in atomic units h — c — e — \)Ss, 



E„i^ — {l{l + l)do- 



m 



{Vol 5^) (n + / + l) 



(n + Z + 1) 



, n = 0, 1,2, •••and Z = 0,1,2, 



(67) 



where n and I are vibrational and orbital quantum numbers, respectively. Also, the wave 
functions become 



Rni{r) = anir-^ | exp(-\/ -2ME„, + ^ii±il^r) 



12 



X (1 - exp(-5r))'+' P„V ^ (1 - 2 exp(-5r)) 



1 



-r ^ exp 



2ME,a , K/ + 1) 



12 



r j (1 — exp(— 5r)) 



2MEni , /(/ + 1) 



+ 



12 



+ ^ + 1 ;l + 2 



2MEni , /(/ + 1 



+ 



12 



-;exp(-(5r) | , 
(68)' 



where E^i </(/ + !) (5^/(24M). The traditional approximation {do = 0) gives 



Enl — — : 



2m 



m 



in + 1 + 1) 



{n + l + l) 



and 



Rni{r) = a„^r-^ (exp(-£„,r)) (1 - exp(-(5r))'+^ p^2s„,/5,2;+i)(^ - 2exp(-(5r)) 
= Qn /^^"'^'^."^ ^^" r~^exp (-enir/5) (1 - exp(-5r))'+^ 
X 2i^i (-^, n + 2 {sni/S + l + l)]l + 2eni/S; ex.p{-5r)) , 
where Eni — \/—2MEni, Eni < for bound state solution. 



(69) 



(70) 
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V. DIRAC EQUATION WITH COUPLING TO r'^ SINGULAR ORBITAL TERM 



In the previous section we have found that the physical quantities hke the energy spectrum 
are critically dependent on the behavior of the system near the singularity. That is why, for 
example, the energy spectrum depends strongly on the angular momentum, which results 
from the singularity of the orbital term, even for high excited states. Since the 
orbital term is too singular, then the validity of all such approximations is limited only to 
very few of the lowest energy states. Therefore, to extend accuracy to higher energy states 
one may attempt to utilize the full advantage of the unique features of Dirac equation. For 
example, the advantage of the Dirac equation over the Schrodinger-like equation is that the 
spin-orbit angular momentum singularity is which is less singular than r~^. Therefore, 
it is more fruitful to perform the analytic approximation of the orbital term in the Dirac 
equation itself, which is a first-order differential equation, not in the resulting second-order 
differential equation. The advantage is that in such case the orbital term is less singular 
since it goes like not like r~^. Therefore, one would expect that the solution of the 
Dirac equation is more accurate by approximating the less singular distribution r~^, which 
makes it possible to extend the validity of the results to higher excitation levels giving better 
analytic approximation for a wider energy spectrum [38] . 

Approximating the spin-orbit term in the Dirac equation (Eqs.(18) and (19)) by a 
function, say, W{r) ^ results in the following second order differential equations that 
should replace Eq. (22) and Eq. (23) (in the relativistic units h—c—1), respectively [38] 

" "'^'^'^ + " ^ ^^"-^ + ^^'^ + ^"'^ " ^0 ^"-^^^^ = °' ^^^^ 

and 

" ^'^'^''^ ~ " ^ ^^"'^ " ^''^ ^^'^ + ^^'^ " ^^'^) ^"'^^'"^ = °' ^^^^ 

where Mg — M — Cg and Mpg — M + Cpg. Note that the resulting proper approximation for 
the term is not as trivial as one would think. That is, the approximation for this term is 
not simply W^{r) but also includes the derivative dW{r)/dr giving the supcrsymmctric form 
W'^{r) ± W'{r). To obtain an alternative solution using this suggested proper approxima- 
tion scheme, we may consider W{r) — 5/ ((e^'' — l) , which is proportional to the Hulthen 
potential. Therefore, we have applied the following proper approximation introduced very 



18 



recently by Alhaidari [38] 



= K^W^r) T ^iW'ir). (73) 



A. Spin Symmetry solution 

We start by solving the Dirac-Hulthen problem in the presence of spin symmetry. We 
approximate the orbital term by a singular function W{r) under certain approximation 
condition that will be maintained throughout the subsection. If we define the variable 
X = e~^^ G [0,+l] and inserting V{r) — —Voe~^^/{l — e~^^), then the positive energy 
Schrodinger-like equation (71) in the new variable x reads as follows 

. l-x d -{(3l + al + K^)x^+{/3l + 2al-K)x-al \ 
dx^^ x{l-x)dx^ a;2(l-a;)2 jPn.[x)-0, [(^) 

provided that 

ai = ^y/M^-E^^, (3, = ^V2(K« + M,)yo. (75) 

Therefore, real solutions are possible only for {Enul < Ms and potential strength Vq > 
(i.e., bound states). Following the procedures explained in the previous section, we can find 
the parametric constants of the NU as fisted in Table 3. Further, the energy equation can 
be obtained with the help of Table 3 and Ref. [22] as 

^/(M - - El^ = ^ [2 (£;„. + M,) V, + kH'] - ^, (76) 

where N2 is given in (57b). The above energy equation has the following simple energy 
spectrum formula 

Kk- TTFT^ \hrF^ + -FT^ y'') 



with 



Qs = H [5'(A^' - «') - 2VoMs] , (78a) 

Ps = V^ + 5^Nl (78b) 

Ws = Ms {MsPs + Qs) + W (2^' - >^^) - ^t] , (78c) 

, n + Z + 1 for K > , , 

Ns='{ , n = 0,l,2,3,--- , (78d) 
n — I for K < 
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where + 4:PsWs > for real spectrum (bound states). For numerical work, Eq. (77) and 
Eq. (78) are used to calculate a few positive energy levels for various values of quantum 
numbers n and I in the spin symmetry limit. We present some numerical values in Table 
2 taking the following values of parameters: M = 5.0 fm^^, Vq = Sq = 3.40 fm^^ and 
Cs = 4.90 fm~^ for the sake of comparison with the previous approximation results 
in Section IV. In referring to Table 2, it should be noted that the spectrum, in the 
approximation scheme, is wide with a fast transition toward the positive energy sector, 
however, in the case of approximation scheme, it is narrow with a slow transition toward 
the positive energy. 

Next, we calculate the upper component wave functions of Eq. (71) in the form of 
hypergemetric function. Moreover, the nonrelativistic limit, is obtained from Eq. (71) by 
setting Ms ^ M {Cs ^ 0), K ^ I > 0, En^ + M ^ 2m, En^-M ^ Enu 2\^ V. Therefore, 
Eq. (76) becomes 



Er,l — 



2 -m{VQ/5^) +1^ (n + Z + l)'" 



(n + / + l) 2 ' ^^^^ 

which is identical to Eq. (14) of Ref. [38] found for V{r) = Vq/ (e'^'" — l) . For the 5*- wave 
{I — 0) restriction of (79) reproduces the well-known nonrelativistic exact result [17]. The 
P term is completely missing from the spectrum formula (67) because of the approximation 
used for and being substituted by another term 5'^l{l + l)/(24m). 
Firstly, the weight function [33] reads 

p(x) =a;2"i(l-a;)2''+\ (80) 

which gives the first piece 

yn.{x) = p(2"i.2.+i) (1 _ ^ ^F, (-n, n + 2 (ai + + 1) ; 1 + 2a,; x) , 

(81) 

and further the second piece reads 

$(x) = (1 - x)^ , 7 = K + 1 (82) 

where a2 and 7 are real positive parameters. Finally, using Eq. (2), we can combine the 
two pieces as 

Fnn{r) = A„,e-^"^^ (1 - e-'^'-)"'"' ^F, (-n, n + 2 (ai + « + 1) ; 1 + 2a,; e"*^) , n = 0, 1, 2, • • 

(83) 
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where is the normahzation factor and k — I ior k > and k — —{l + l) ior k < 0. 
The bound state solution requires that the hypergeometric series terminate. The lower 
component wave functions are calculated from Eq. (18) as 

= K:Tm: \ il-e-sr) - + -r) 

n5 (n + 2q;i + + 2) _ 5 +1 . ai+i 
^^"'^ (K. + M,) (1 + 2a0 ^ ^ i J 

X 2F1 (-n + l;n + 2(Q;i + K + l) + l;2(l + Q;i);e-'''') . (84) 
Moreover, from Eq. (83), the nonrelativistic radial wave function reads 

Rni{r) = Ar,ie-'^-^'{l-e-'')'^' ^F^ (-n, n + 2 (A^^ + / + 1) ; 1 + 2An;; e"-^^) , n = 0,l,2,-- 

(85) 

where Xni — y/ —2'mEni/ 5, E^i < which is defined via (79). 
B. Pseudospin symmetry solution 

The pseudospin symmetry solutions could simply be found by applying the following map 
on the spin symmetry solution in the previous section (both the energy spectrum and the 
spinor wave functions) [35,38] 

Fnn{r) O GnK{r), Enn ^ -^n«, « -K, Vq ^ -Vq, Cs ^ Cps, (86) 

which leads to the generation of Eq. (72) from Eq. (71). Making the change of variables, 
X = e~^'^ e [0, +1], we can rewrite Eq. (72) as 

+ + 2n ^2 -]Gr^,{x) = 0, 87 

ax^ x[l — X) ax a;^(l — xj"^ 



provided that 



Note that the parametric constants for the present case are hsted in Table 3. The energy 
equation becomes 

^(M+^TT^ = ^ [2 {E^, - M - C,,) Vo + k'S'] - (89) 



21 



where A^i is given in (35). The above energy equation has the following simple energy 
spectrum formula 

p± _ Qps , / Qps , , , 

Kk - 775- ± \ + -5-' (y") 



2P A/ 4P2 P ' 



with 



Qps = Vo [6\Nl, - K^) + 2VoMps] , (91a) 
Pps^V,' + 5'N^,, (91b) 
VTp, = Mp, {MpsPps - Qps) + [/^' i2N^s - - Ks] ' (91c) 

{n + r+ 1 for K>0 
, n = 0,l,2,3,--- , (91d) 
n — i for K < 

where + 4PpslVp5 > for real spectrum (bound states). Tables IV and V give approxi- 
mation to the exact spin and pseudospin symmetry, respectively. Further, we have provided 
two different approximation models ( i.e., and r"^) for the sake of comparison. The 
approximation scheme (either conventional or improved) is found to be more sensitive to 
spin-orbit quantum number k than the r"^ proper approximation model [38]. The latter is 
found to be valid for wide range energy spectrum (see Tables II, IV and V).[l] 
The lower spinor can be found as 

Gnn{r) = Bn^e-'"'''- (l - e-^')" ^F, (-n, n + 2 {a2 + k) ;! + 2^2; e"-^^) , n = 0, 1, 2, ■ • • . 

(92) 

where is the normalization factor, k — l + lfoTK>0 and k — —I for k < 0. The upper 
spinor component wave functions are calculated from Eq. (19) as 

""-^'^ ^ ( (T^ - - ^) ""-^'^ 

,0 n(5 (n + 2a2 + 2k) /,-^r^"2+l 
+''"'^(M,,-E..)(l + 2a2)^ ^ j 

X 2F1 (-n + l;n + 2(Q;2 + K) + l;2(l + Q;2);e-^''). (93) 



[1] We have found a considerable discrepancy in the numerical results of the two approximation schemes. 
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VI. RESULTS AND CONCLUSIONS 



In the pseudospin symmetry case, the variation of the positive (negative) energy spectrum 
(£'"-) with the screening parameter 5 is shown in Fig. 2a (Fig. 2b), with a suitable set 

Tl/t Tit 

of physical parameter values Cps = —4.90 /m~^, M — 5.0 fm~^ and Aq = 3.40 fm^^. For 
specific values of quantum numbers n and it is noted that when the screening parameter 
S increases, the positive energy increases in the negative energy direction ~ 0.03 fm~^, 
~ 0.13 /m~^ and ~ 0.35 fm~^ for pseudoorbital quantum numbers I — 1, I — 3 and / — 5, 
respectively, with a small energy difference between states is small (~ 0.02—0.04 fm~^) when 
5 — 0.20. However, the negative energy spectrum increases in the negative energy direction 
~ 0.50 fm~^, ~ 1.40 fm~^ and ~ 1.90 fm~^ lor I — 1, I — ?, and I — 5, respectively, with 
energy difference (~ 0.45 — 0.50 fm~^) when 5 = 0.20. For example, when 5 = 0.20, we have 
~ 0.01, 0.03 and 0.08 fm'^ and E^j ~ 0.15, 0.60 and 1.15 fm^^ with 7'= 1,7'= 3 and 
I = 5, respectively. This large spacing returns to the new shifting energy term 5^/(Z + 1)/12. 
A more strongly binding energy occurs for E^^ when < 5 < 0.1 (Fig. 2a) but for when 
the screening parameter is lower, i.e., < 5 < 0.05 (Fig. 2b). 

In the spin symmetry case, the variation of the energy spectra {E^^ and £'~^) with the 
screening parameter 5 is shown in Figs. 2a and b, with a suitable choice of physical parameter 
values Cs = 4.90 fm'^, M = 5.0 fm~^ and Sq = 3.40 fm~^. The positive (negative) energy 
spectrum is plotted in Fig. 3a (Fig. 3b). For specific values of quantum numbers n and «;(/), 
it is noted that when the screening parameter 5 increases, the positive energy increases with 
a large amount in the positive energy direction (~ 4.6 fm~^) and difference between states 
is large (~ 1.3 — 2.0 fm~^) whereas the negative energy spectrum increases with small 
amount in the positive energy direction (~ —0.10 — 0.05 fm~^) and difference in energy 
spacing is (~ 0.55 — 0.85 fm~^). For example, when 6 = 0.20, we have E^^ ~ 0.2, 1.2 and 

2.4 /m-i and E"^ 0.095, -0.085 and -0.065 fm'^ for the orbital states / = 1, 3 and 5, 

respectively. These large (small) shifts return to the new shifting energy term 5^l{l + 1)/12. 
A more strongly binding energy occurs for when < S < 0.05 (Fig. 3a) but for £^~^ 
when the screening parameter is higher, i.e., < 5 < 0.1 (Fig. 3b). 

Under the pseudospin symmetry, the energy-mass curves are plotted versus mass for 
the pseudospin orbital quantum numbers / = 1, Z = 3 and / = 5 by taking the pseudospin 
constant Cpg = —4.90 fm~^ and the screening parameter 5 = 0.25 for a given radial quantum 
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number n — as shown in Fig. 4. There are two different regions of energy spectrum (E'^^ 
and E~i^) with the mass as shown in Fig. 4. In the positive energy part, nearly in 
the region < M < 2.4 fm~^, the energy spectrum is all in the negative region and the 
energy increases in the direction of the positive energy as / increases. In the negative energy 
part, E~^, nearly in the region 4.3 < M < 12 fm"^, the energy spectrum decreases in the 
direction of the negative energy when I increases. Furthermore, under the spin symmetry, 
the energy-mass curves are plotted versus mass for the orbital quantum numbers 1 — 1,1 — 3 
and I — 5hy taking the spin constant Cg — 4.90 fm~^ and the screening parameter 5 — 0.25 
for a given radial quantum number n = as shown in Fig. 5. There are two different regions 
of energy spectrum {E^^ and E^^) versus mass as shown in Fig. 5. In the negative energy 
part, E~^, nearly in the region < M < 0.6 fm^^, the energy spectrum is in the positive 
region. The energy decreases in the direction of the negative energy as I increases. In the 
positive energy part, E^^, nearly in the region 3.0 < M < 12 fm~^, the energy spectrum 
increases in the direction of the positive energy when I increases. 

In Fig. 6, we have plotted the energy spectrum versus the pseudospin constant Cps 
for the parameters values M = 5.0 fm~^, Aq = 3.40 fm~^ and S = 0.25. The negative 
values of Cps show more strongly binding energies for Cps < —10 fm'^ in the and less 
strongly binding energies for Cps > — 6 fm~^ in the for all I values. The energy for the 
constants —9 fm~^ < Cps < —5 fm~^ still show the negative energy up to the zero axis. 
Furthermore, in Fig. 7, we have also plotted the energy spectrum versus the spin constant 
Cs for the parameters values M — 5.0 fm~^, Eq = 3.40 fm~^ and 5 — 0.25. The positive 
values of Cs show more strongly binding energies for Cs > 15 /m~^ in the E~^ and less 
strongly binding energies for Cs < 10 fm~^ in the E^^^ for all / values. The energy for the 
constants —2 fm~^ < Cs < —20 fm^^ still show the negative energy up to the zero axis. 
For the case considered in Fig. 2a (Fig. 2b) where Cps — —4.90 fm~^, M — 5.0 fm~^ and 
Aq = 3.40 fm"^ in the less strongly binding energies shows the negative energy up to the 
zero axis and falls in the region However, the case where Cg — 4.90 fm~^, M — 5.0 
fm~^ and Sq = 3.40 fm"^ considered in Fig. 3a (Fig. 3b) falls in the less strongly binding 
energies in the region E^^. 

We have seen that the Dirac equation for the Hulthen potential based on spin symmetry 
and pseudospin symmetry limitations can be solved approximately for any arbitrary spin- 
orbital K state within the framework of the Dirac theory. By using the basic ideas of 
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the parametric generalization of the NU method, the approximated positive and negative 
energy eigenvalues for the arbitrary spin-orbital (pseudo-orbital) angular momentum I {I 
) are obtained. An improved approximation scheme is used to deal with the centrifugal 
/(/ + l)/r^ (pseudocentrifugal /(/ + 1)/t^) potential term. The energy spectrum for any / (/) 
states is obtained analytically under the spin symmetry, A(r) = (pseudospin symmetry, 
E(r) = 0) limitations, the energy relations in the Dirac equation with equal scalar and vector 
Hulthen potentials are recovered to see degenerate states. The relativistic bound state energy 
eigenvalues and the correspondinf two-component spinor wave functions have been easily 
reduced to the non-relativistic limits by applying appropriate parameters replacements. 

Finally, it is noted, from Tables 2, 4 and 5, that analytic solution of the Dirac equation is 
more accurate by approximating the less singular distribution spin-orbit angular momentum 
term r~^, which makes it possible to extend the validity of results to higher excitation 
levels giving better analytic (numerical) approximation for a wider range spectrum since 
the dependence of r"^ on the angular quantum number is less than the too singular term r"^. 
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FIG. 1: A plot of the variation of the singular orbital term (dotted-solid line) with (a) the 

approximation of Ref. 4 (dash line) , the conventional Greene- Aldrich of Ref . 36 (dash-dot line) and 
improved (solid line) approximations of with respect to r where S = 0.1 fm'~^, (b) improved 
approximation with various shifting constants. 



FIG. 2: The pseudospin symmetry energy spectrum versus the screening parameter 6 for (a) 
and (b) E~i^ in units of fm~^. 



FIG. 4: The pseudospin symmetry energy spectrum versus the mass. The solid, dash and dash- 
dot lines of the energy spectra including £'+^ and E~i^ are plotted for Z = 1, Z = 3 and 1 = 5, 
respectively. 



FIG. 5: The spin symmetry energy spectrum versus the mass. The solid, dash and dash-dot lines 
energy spectra including and E~i^ are plotted for Z = 1, Z = 3 and I = 5, respectively. 



FIG. 3: The spin symmetry energy spectrum versus the screening parameter 6 for (a) E^^ and (b) 
E~ in units of fm^^. 
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FIG. 6: The pseudospin symmetry energy spectrum versus the constant Cps- The sohd, dash and 
dash-dot Unes energy spectra including and are plotted loi I = 1, I = 3 and I = 5, 
respectively. 



FIG. 7: The spin symmetry energy spectrum versus the constant C^. The solid, dash and dash-dot 

lines energy spectra including E^^ and are plotted for / = 1, / = 3 and / = 5, respectively. 



TABLE I: The negative bound state energy eigenvalues in units of fm ^ of the pseudospin sym- 
metry Hulthen potential for various values of n, Z = Z -|- 1 and S. 



I n,K < 0,K > Q S States Enn Enn [4] I n,K<0,K>Od States E, 



1 1, 


-1,2 


0.025 


(1^1/2,04/2) 


0.0972235 


0.0963638 1 


2, 


-1,2 


0.025 


(251/2,14/2) 


0.0938034 






0.100 




0.0561798 


0.0425738 






0.100 




0.0038600 






0.175 




-0.0302923 


-0.0710009 






0.175 




-0.175897( 






0.250 




-0.1544010 


-0.2346580 






0.250 




-0.412557( 


2 1, 


-2,3 


0.025 


(1^3/2,0/5/2) 


0.0937343 


0.0912282 2 


2, 


-2,3 


0.025 


(2^3/2,1/5/2) 


0.0889591 






0.100 




0.00275013 


-0.0363590 






0.100 




-0.0673921 






0.175 




-0.1793260 


-0.2930130 






0.175 




-0.359049( 






0.250 




-0.4196540 


-0.6351320 






0.250 




-0.704102( 


3 1, 


-3,4 


0.025 


(1^^5/2,057/2) 


0.0888560 


0.0839128 3 


2, 


-3,4 


0.025 


(24/2,157/2) 


0.0827390 






0.100 




-0.0690512 


-0.1447100 






0.100 




-0.154261( 






0.175 




-0.3642070 


-0.5760950 






0.175 




-0.561156( 






0.250 




-0.7148860 


-1.0984500 






0.250 




-0.987242( 


4 1, 


-4,5 


0.025 


(1/7/2,0/19/2) 


0.08260190 


0.0744360 4 


2, 


-4,5 


0.025 


(2/7/2, 1^9/2) 


0.0751593 






0.100 




-0.1564720 


-0.2784550 






0.100 




-0.253646( 






0.175 




-0.5680850 


-0.8953110 






0.175 




-0.767387( 






0.250 




-1.0019200 


-1.5671200 






0.250 




-1.238430( 
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TABLE II: The positive bound state energy eigenvalues in units of fm of the spin-symmetry 
Hulthen potential for various values of n, I and 6. 



I n,K<0,K>0 S States (present)'' (present)*/ n,K<0,K>0 S States Ej, 



1 n 


— 9 1 




1 1 I'M- j~ 1 ln„ i~ 1 

(,UPl/2) ^P3/2) 


— n nQ49nn^ 


— n nQQ'^Q1 


1 
1 


1 


— 9 1 


n n9'i Ht), /„ in„ /^^i 

U.UZO \i.Pl/2i '-P3/2) 


— n 

— u. 






inn 




— n nnszino^'i 


— n nQ^'in9'i 








n inn 

U. -LUU 


-un 






0.175 




+0.1727090 


-0.0803626 








0.175 


+0. 










-un 4^'^fi'^nn 


— n nfin74/i7 








n 9'in 

u.zuu 


-un 


9 n 


— ^ 9 












z 


1 


— ^ 9 




— n 

— u. 






n inn 




-un in97fi^n 


— n ny^nynzi 








n inn 

u. ±uu 


4-n 






0.175 




+0 4840740 


-0 0249639 








0.175 


+0 






n 9^n 




^n QQ1 "^(^Rn 


-un nziQi (^n^ 








n 9^n 

u.zou 




^ n 




n n9^ 




— n nvfis^ns 

— u.u / uoouo 


— n nnvnAQi 


Q 


1 

J-, 




n n9f^ n ,„ 1 f» 

U.UZO ;^i/5/2, ^J7/2) 


— n 

— u. 






n inn 




_i_n 










n inn 

u. ±uu 


-un 






0.175 




+0.8720970 


+0.0385481 








0.175 


+1. 






0.250 




+1.6200500 


+0.1706690 








0.250 


+2. 


4 0, 


-5,4 


0.025 


(057/2,059/2) 


-0.0638592 


-0.0952974 


4 


1, 


-5,4 


0.025 (157/2,159/2) 


-0. 






0.100 




+0.4345750 


-0.0262998 








0.100 


+0. 






0.175 




+1.3032300 


+0.1159560 








0.175 


+1. 






0.250 




+2.2434000 


+0.3130470 








0.250 


+2. 



"Improved approximation to a more singular orbital term r 
''Proper approximation to a less singular term [38]. 
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TABLE III: Specific values of the NU constants based on the spin and pseudospin symmetric Dirac- 
Huhhen problem considering the recently introduced proper approximation to the less singularity 
orbital term. 



Spin symmetry: 


Pseudospin symmetry: 


Cl = 1 


Cl = 1 


<-o — 1 

C2 — i 


C2 — i 


^ 1 
C3 = 1 


^ 1 
C3 = 1 


C4 = 1 


C4 = 1 


C5 = 


C5 = 


C6 - -2 


C6 - -2 






^ „ 2 /o2 1 

C8 = — 2a;^ — Pi + K 


„ 0^2 02 


Cg = af 


eg = 02 


cio = 3(2«; + 1) 


1 1 \2 

Clo = 2(2k - 1) 


cii = 2ai 


Cii = 2q;2 


C12 = 2k + 1 


C12 = 2k — 1 


ci3 = ai 


Cl3 = a2 


Ci4 = At + 1 


Ci4 = K 


Ci5 = 2k + 1 


C15 = 2k — 1 


C16 = + 1 


CI6 = K 


^1 = af + /3? + 


^1 = q2 + /3i + k2 


^2 = 2a2 + - K 


^2 = 2a2 + /32 + K 


^3 = a? = <5-^ (M| - 


EL) ^s = al = 6-^M^,-. 
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TABLE IV: Approximation of the negative bound state energy eigenvalues based on the exact 
pseudospin symmetry {Cps = 0) Hulthen potential for various values of n, I = i + 1 and S. 



I n,K<0,K>0 6 States (present)'* (present)* I n,K < 0, k > Q 6 States E^ 



1 1 

-"-J 


— 1 9 




\^^1/2j ^"'3/2) 


4 QMOS 


4 QQfil 1 


1 9 


— 1 9 


U.UZO (^Z62/2) ^"'3/2) 


4 Q7 






inn 






4 Q^891 






n inn 








0.175 




4.28511 


4.81377 






0.175 


3.81 












4 fi9Qnfi 






n 9'in 


9 SCI 


9 1 




n 09^ 


l-'^P3/2) 


4 Q71 fi^ 


4 QQ97n 


9 9 


—9 ^ 
z, 


n n95 Oti^ 1 fr- 

U.uzo (,zp3/2, ±/5/2j 


4 






n inn 




4 ^fiSS^ 
^.ouooo 


4 884fiQ 






n inn 


4 ^4 






0.175 




3.80980 


4.65663 






0.175 


3.28 






n 9'in 




9 SQ'^ni 


4 '^97Q9 






n 9'in 








n n9^ 


\^'^b/2->^97/2) 




4 Q8891 


^ 9 
z. 


— ^ 4 


n n9'i (9ri^ 1 /^^i 

u.uzo (,za5/2, iy7/2j 


4 






n inn 




4 3455fi 


4 81 51 5 






n inn 


4 OQ 






0.175 




3.28126 


4.45771 






0.175 


2.73 






0.250 




2.12740 


3.96084 






0.250 


1.42 


4 1, 


-4,5 


0.025 


(1/7/2,0/19/2) 


4.93644 


4.98265 


4 2, 


-4,5 


0.025 (2/7/2,1/19/2) 


4.91 






0.100 




4.08954 


4.73030 






0.100 


3.80 






0.175 




2.73540 


4.22266 






0.175 


2.20 






0.250 




1.42282 


3.54589 






0.250 


0.81 



"Improved approximation to a more singular orbital term r 
''Proper approximation to a less singular term [38]. 
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TABLE V: Approximation of the positive bound state energy eigenvalues based on the exact spin 
symmetry (C^ = 0) Hulthen potential for various values of n, / and S. 



I n,K<0,K>0 S States E^^k (present)'* (present)^/ n,K<0,K>0 5 States E„ 



1 





—2 1 


0.025 


\yyi/2i 


—4 98993 


—4.99731 


1 1 


—2 1 


025 f1ni /o Ino/nl 
u.u^u (,-i-/-'i/2) ^i'Z/2) 


—4. 








ino 






—4 9571 8 






100 


—4. 








0.175 




-4.52642 


-4.86979 






0.175 


-3. 








2.50 




—4 07294 


—4 7371 7 






250 




2 


n 


— S 2 


025 


V""'3/2) 


—4 97737 


—4 99356 


2 1 


—3 2 


025 (^(Jr,ln Ir/c/o) 


—4. 








100 




—4 fi481 5 


—4 89773 






100 


—4. 








0.175 




-3.98590 


-4.69175 






0.175 


-3. 








2,50 




— 3 1 031 7 


—4 38588 






250 


— 2. 


3 


n 


—4 3 


0.025 




—4 95982 


—4.98847 


3 1 


—4 3 


025 fl /"c /n 1 /■^/n'l 


—4. 








0.100 




—4 38880 


—4.81796 






0.100 


—4. 








0.175 




-3.31174 


-4.45782 






0.175 


-2. 








0.250 




-2.00840 


-3.93859 






0.250 


-0. 


4 


0, 


-5,4 


0.025 


(057/2,05(9/2) 


-4.93733 


-4.98196 


4 1, 


-5,4 


0.025 (157/2,159/2) 


-4. 








0.100 




-4.07241 


-4.71713 






0.100 


-3. 








0.175 




-2.56164 


-4.17002 






0.175 


-1. 








0.250 




-0.91879 


-3.40926 






0.250 


+0. 



"Improved approximation to a more singular orbital term r 
''Proper approximation to a less singular term [38]. 
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